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Geometrically frustrated materials have a ground-state degeneracy 
that may be lifted by subtle effects, such as higher order inter- 
actions causing small energetic preferences for ordered structures. 
Alternatively, ordering may result from entropic differences between 
configurations in an effect termed order-by-disorder. Motivated by 
recent experiments in a frustrated colloidal system in which ordering 
is suspected to result from entropy, we consider in this paper the an- 
tiferromagnetic Ising model on a deformable triangular lattice. We 
calculate the displacements exactly at the microscopic level, and 
contrary to previous studies, find a partially disordered ground state 
of randomly zigzagging stripes. Each such configuration is deformed 
differently and thus has a unique phonon spectrum with distinct en- 
tropy, lifting the degeneracy at finite temperature. Nonetheless, due 
to the free-energy barriers between the ground-state configurations, 
the system falls into a disordered glassy state. 
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Frustrated systems are characterized by interactions that 
may not be satisfied simultaneously pQ. This leads to 
a degenerate and thus disordered ground state, and naively 
one would expect to have disorder down to zero temper- 
ature [2]. However, the frustrated phase is very sensitive 
to small perturbations that can order it. These include 
anisotropic [3] or longer-range [4] interactions and lattice de- 
formations 5, 6, 7 8, 9, 10 as well as entropic effects that may 
lift the ground-state degeneracy at finite temperature, in a 
process termed order-by-disorder [11M12I [T51 Tl4 15, 16]. Frus- 
tration and its relief due to order-by-disorder are traditionally 
investigated in antiferromagnets and, in particular, in com- 
pounds that have a triangular lattice structure [TTl fl8l [T9l [20] . 
Recent experiments have demonstrated that artificial systems 
made of mesoscopic building blocks such as single-domain 
magnetic islands [21] or colloidal spheres [22] exhibit behav- 
ior which is similar to that of magnetic systems comprised of 
atomic-scale particles. Such mesoscopic systems enable direct 
visualization of the dynamics at the single-particle level and 
thus provide insight into the microscopic physical mechanisms 
responsible for the peculiar properties of frustrated matter. In 
this paper we address the theoretical underpinnings of these 
new experimental systems and thus provide a vital step in 
understanding their unusual behavior and its connection to 
atomic-scale antiferromagnets. 

A densely-packed monolayer of hard spheres buckles out 
of its plane when it is confined between walls that are sep- 
arated by slightly more than a single sphere diameter [23] . 
Entropic forces depending only on geometry give rise to 
effective anti-ferromagnetic interactions favoring motion of 
neighboring spheres toward opposite walls |24| . leading at 
high densities to stripes of alternating up and down spheres. 
The close-packed state is highly degenerate: the same max- 
imal density is obtained by straight stripes or by any set 
of parallel stripes that zigzag within the hexagonally-packed 
layer [22] [24] [55] [26] [27] [28] . Recent experiments indicate a 



possible preference of the stripes to be straight rather than 
to zigzag randomly in the plane, suggesting that at densities 
below close-packing there is an order-by-disorder effect giving 
an entropic advantage for straight stripes [22] . 

This quasi-two-dimensional problem is strongly related to 
the old, yet unsolved question of what is the stable high- 
density structure of hard spheres, face-centered-cubic (FCC) 
or hexagonal-close-packing (HCP) 29, 30, 31 . In three di- 
mensions, maximal density is obtained by stacking hexago- 
nally packed layers with arbitrary sideways shifts between the 
close-packed positions. And, as for the buckled monolayer, 
slightly below close-packing it is not clear which structure 
has the greatest entropy. Experiments on colloidal crystals 
that were grown slowly enough exhibit FCC order [32], while 
there is controversy on whether the theoretical estimates have 
reached the accuracy required to resolve the elusive entropic 
difference between FCC and HCP [331 [34] 1351 l36l [371 [38] . 

As for the FCC vs HCP question, the entropy, or free vol- 
ume, of buckled hard spheres is a collective function of the 
positions of all particles, thus the ability to obtain analyti- 
cal results for it is very limited. Thus, instead of approxi- 
mating the entropy of this hard-sphere system, we consider 
the anti-ferromagentic Ising model on a deformable triangu- 
lar lattice. This model has the same degenerate ground state 
of zigzagging stripes as the colloidal system, and for it we 
can exactly calculate the free-energy difference between the 
competing configurations. We find that straight stripes are 
always favored entropically. However, the free-energy barriers 
between various ground states are huge compared to this en- 
tropic gain, causing the system to fall into a disordered glassy 
state upon cooling. 

In our Ising model, each site i at continuous position r; on 
a triangular network is occupied by a particle of discrete spin 
<Ji = ±1, and each nearest neighbor bond comprises a har- 
monic spring of stiffness K and relaxed length a, see Fig. [TJ\. 
The internal energy of the system depends only on the spin 
product OiOj of neighboring particles and on their relative po- 
sitions, 8r\j = \ fi — fj\ — a. The model's Hamiltonian is given 
by the following sum over all nearest-neighbor pairs (ij), 



(1 — eSnj) JoiOj + —5r% 
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The antiferromagnetic interaction strength is equal to J > 
for particles separated by the relaxed spring length a, and 
decreases linearly with distance at a rate e > 0. This mim- 
ics the buckled colloidal monolayer in the following way: For 
hard spheres, free energy is determined by their free volume. 
However, it is convenient to reduce this to an effective two- 
particle repulsive interaction. While the particles are confined 
to be quasi-two-dimensional, they do have a limited freedom 
to move in the vertical direction, and thus nearest neighbors 
prefer to sit at opposite heights. We can now map the verti- 
cal position of each particle onto an Ising degree of freedom, 
due to a preference for either an "up" or "down" state. This 
leads to an effective two-dimensional antiferromagnetic Ising 
model, with the coupling between the elastic and the mag- 
netic degrees of freedom as reflected in our Hamiltonian {l}, 
i.e. the effective antiferromagnetic interaction decays as the 
in-plane separation increases [24] , 

We model the in-plane entropic repulsion between neigh- 
boring spheres by a rotationally invariant central-force har- 
monic spring [Eq. ([TJ] . Approximate versions of this model 
have been studied previously: Ref. [6] uses a linearized long- 
wavelength elastic energy, invariant with respect to only in- 
finitesimal rotations, that limits deformations to be small and 
to vary slowly in space, while in Ref. [5] deformations are 
taken to be uniform. 

In this paper we present results for the case in which the 
total area of the system is fixed to the area of a triangular 
lattice with lattice constant equal to a. However, we obtained 
qualitatively similar results for systems compressed or dilated 
with respect to this simple situation. From the equivalence be- 
tween ensembles we thus expect to get the same results also 
when considering the fixed pressure case, which is probably 
more appropriate for the colloidal experiments. 

As in the rigid triangular-lattice model 2,3, the anti- 
ferromagnetic interactions along each three-particle loop in 
our deformable network cannot be satisfied simultaneously, 
and energy is minimized by having two satisfied (oiOj = — 1) 
bonds and a single frustrated (oiOj = 1) bond around each 
triangular plaquette. Because of the magnetoelastic coupling, 
energy may be lowered by stretching the frustrated bonds by 
a factor / and compressing the satisfied ones by a factor s. 
We fix the area of each resulting isosceles triangle to be that of 
the initial equilateral plaquette with sides a, and thus we can 
parametrize the deformation by the head angle /3 (see Fig. [1)3, 
inset): /(/?) = 3 1/4 (tan f ) 1/2 , s(/3) = 3 1/4 (2sin/3)- 1/2 . Min- 
imizing the energy fl]) with respect to /? yields (see SI Ap- 
pendix) 

(2s' - f')Je + [2(s - l)s' + (/ - 1)/'] Ka = 0. [2] 

Figure [TJ3 shows how the triangles deform from /3 = 60° to- 
ward /3 — 180° as the ratio b = -j^ of the magnetoelastic 
interaction strength to the lattice rigidity grows. 

Thus, each plaquette of the triangular lattice would min- 
imize its energy by deforming into an isosceles triangle. We 
now show how this can be accommodated in the ground state 
by global deformations of the system. In the rigid triangular 
lattice, since each triangle must have exactly one frustrated 
bond, five nearest-neighbor configurations are allowed in the 
ground state (Fig[Tp) [39| . Requiring that the angle oppos- 
ing each frustrated bond deforms to /3 > 60° selects config- 
urations (Hi) and (iv) |22l 124] . which give rise to zigzagging 
stripes. Thus, selecting a stripe of frustrated bonds defines a 
ground state of the lattice. It can be constructed by starting 
with a row of alternating spins (along the horizontal axis) and 
stacking copies of this row (along the vertical axis), as shown 
in Fig. [2] The intra-row and inter-row separations are sa and 



sasin/3, respectively, and the lateral shifts are determined by 
the arbitrary polarity of each row with respect to the one pre- 
ceding it. Since these rows of alternating spins may be in any 
of the three principal directions of the network, and for a sys- 
tem of N particles, each of the y~N rows may be in one of two 
states, the ground-state degeneracy is 3 • 2^ '. 

This partially-disordered ground state is realized by defor- 
mations that may vary rapidly in space; if deformations are 
assumed to be homogeneous [5] or to vary slowly [BJ, straight 
stripes are selected. In particular, the zigzagging stripes that 
minimize our microscopic Hamiltonian ([1} have a higher en- 
ergy than straight stripes in the coarse-grained Hamiltonian 
considered in [6]. The ground state of zigzagging stripes that 
we find here is precisely the state that maximizes the pack- 
ing density of buckled spheres [221 124] . which constitutes the 
connection of our model to the colloidal system. 

Given the high ground-state degeneracy and the fact that 
it takes a discrete energy (of order J) to flip a spin in the 
ground state, one might naively expect to find a stable phase 
of randomly zigzagging stripes at sufficiently low temperature 
(ksT <C J). However, at positive temperature, the entropy of 
particle fluctuations around the energy-minimizing position 
is different in the different ground-state configurations. We 
will show that the state with straight stripes has lower free 
energy (or greater entropy) than other states, making it the 
stable thermodynamic phase at arbitrarily low temperature. 
Nonetheless, the resulting free-energy differences between the 
different zigzagging realizations are much smaller than the 
free-energy barriers between the ground-state configurations 
and the system typically falls into disordered meta-stable con- 
figurations. 

To demonstrate the peculiar slow dynamics of this model, 
we used Monte-Carlo simulations in which individual particles 
can flip their spin or move continuously in the plane and in 
which the simulation box may change its shape to accommo- 
date the global deformations of the network [24] . We start 
from a disordered state at high temperature, and we slowly 
cool the system at a rate r (time is measured in attempted 
Monte-Carlo steps per particle). Figure [3]4. shows that at 
high temperature the system follows an equilibrium curve ir- 
respective of cooling rate, while below a certain temperature 
(T ~ 0.6 for the parameter values shown here), the system's 
energy has a clear cooling rate dependence. Apparently, if the 
cooling is slow enough (r < 10~ 6 here), the system manages 
to reach the ground state. 

Figure [4j\ shows that indeed if the system is cooled too 
rapidly (r = 10~ 4 ) it falls into a disordered state with mul- 
tiple small domains of a local stripy structure. For a slower 
cooling rate (r = 10~ 6 , Fig. [3)3), the system finds a ground- 
state configuration with zigzagging stripes, such that the local 
configurations (in) and (iv) defined in Fig. QJJ are roughly 
equally represented, thus corroborating this phase as ran- 
domly zigzagging stripes. Figure 2P shows that for even 
slower cooling (r = 10 -8 ), there is preference for ground-state 
configurations in which the stripes are more straight than 
bent, namely configuration (Hi) is preferred over (iv). This 
is quantified in Fig. [3] B, C, which shows results obtained by 
averaging over multiple realizations for each system size and 
cooling rate. For very fast cooling, the system remains quite 
disordered with P(iii) / P(iv) ss 0.5, in accordance with the 
combinatoric weights of these two local configurations. For 
r < 10 -6 , the system manages to find its ground state, as 
seen by the fact that P(iii) + P(iv) = 1, namely all particles 
are in one of these two states. For r = 10~ 6 these two states 
are equally probable, P(iii) ~ P(iv), as one would expect for 
randomly zigzagging stripes. However, as the cooling rate is 
decreased even more, there is a clear preference for the straight 
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configuration (Hi) over the bent one (it)) - This is in qualitative 
agreement with recent experiments 40 exploring the cooling 
rate dependence of our colloidal antiferromagnet [22] . 

Although the fsing degrees of freedom in our model do 
not have natural physical dynamics associated with them, the 
Monte-Carlo dynamics we employ are useful in exploring ac- 
tual out-of-equilibrium dynamics of spin systems |41l 1421 1431 
144] . Thus, to our understanding, the cooling rate in our simu- 
lations should be interpreted as being directly proportional to 
the experimental cooling rate, and it would be interesting to 
explore the connection between the transition rates we apply 
in our Monte-Carlo scheme and the actual physical quantities 
in the experiments. 

The preference we observe for straight over bent stripes 
comes from an entropic difference between the two, as is ex- 
plained in the following low-temperature expansion. Expand- 
ing the Hamiltonian about the fixed-spin ground state to low- 
est order, we find terms quadratic in particle displacements. 
Thus, the lowest-energy excitations are harmonic modes of vi- 
brations with frequencies which we denote by {u^} (note that 
k indexes the normal modes and does not necessarily refer 
to a wavevector). We assume the temperature is low enough 
(k B T <C J) to ignore spin flips, yet high enough (k B T 3> Fuo) 
to ignore quantum effects. Therefore, we use the canonical 
partition function of a classical harmonic oscillator of fre- 
quency tj, Z oc k B T/uj, to write the free energy of the system 
as the following sum over all normal modes 



T = -k B T\oe,Z = -k B T 



k B T 



+ const. 



[3] 



We emphasize at this point that the frequencies {uik} refer 
to the fast oscillations of the positional degrees of freedom, 
which at low enough temperature are expected to be much 
faster than the spin flips. 

The different spin states {<Ji} of the various ground-state 
configurations impose different deformations {r*i}, thus each 
one has a distinct spectrum {oJk} of eigen-frequencies, and a 
different entropy, which is expressed as a different tempera- 
ture dependence of the free energy T in Eq. (|3}- At T = 0, 
entropy is irrelevant and all ground-state configurations are 
equivalent, but at T > this degeneracy is lifted and the 
configuration with the minimal free energy is selected. The 
frequencies of all vibrational modes scale as uj'f, oc K/M, where 
M is the mass of each particle. The normalized spectrum de- 
pends on the spin state {ai} and on the deformation angle ft, 
which is set by b = (see Eq. ([2]) and Fig. [1)3). Hence we 
define 



A({oi},p) 




[4] 



and write T = Nk B T [A({<Ji}, /3) + const.]. Thus minimizing 
J- is equivalent to minimizing A. 

We can analytically calculate the normal modes of vibra- 
tions of the deformed lattice for any ground state with a 
periodic repetition of straight (Hi) and bent (iv) segments 
(see SI Appendix). In such cases the normal modes are 
phonons and the index k labeling them may be associated 
with their wavevector. The free-energy coefficient, A, may 
then be obtained by numerically summing over the appro- 
priate Brillouin zone. Figure [2] highlights the unit cells in 
ground states with periods of one, two, and five particles. 
Now we consider ground states with larger unit cells, and 
show that the free energy is mainly determined by the frac- 
tion p s = P(iii) /[P(iii) + P(iv)] of straight segments. To 
clarify which configurations we used in these calculations, we 
show in Fig. [5] several representative unit cells. For example, 



p s — 1/3 in a single configuration with three particles in its 
unit cell (Fig. [5j\), but also in two configurations with six 
particles (Fig. 03), in six with nine particle (Fig. [5p), and in 
an increasing number of other configurations with larger unit 
cells. Similarly, for p s — 1/2, we show the two configurations 
with four particles per unit cell in Fig. [5p and the six with 
eight particles in Fig. [5p. 

Figure [SJ\-C shows the results obtained by numerically 
evaluating Eq. Q for ground states with unit cells consisting 
of up to 13 particles, plotted vs p s (see SI Appendix for more 
details and Tables S4-S5 for the entire data plotted). We find 
that the free energy is bounded between the extreme cases of 
straight stripes (p s = 1, Fig. [2j\) and bent stripes (p s — 0, 
Fig. [2J3)- Moreover, the approximate collapse of the results 
for all calculated unit cells to a single curve provides strong 
support for the hypothesis that only the ratio between straight 
and bent segments is important and not the order in which 
they are positioned or the period of the pattern they form. 
Finally, the roughly linear dependence of A on p s indicates 
that the free energy may be approximated as a linear combi- 
nation of contributions from the straight and bent segments. 
Since the ground state is determined by a one-dimensional 
sequence of straight and bent segments, this result implies 
that the system may be well approximated by an effective 
one-dimensional non-interacting Hamiltonian. For j3 — 155° 
(Fig. HJP) we observe substantial deviations from the predic- 
tion of such a non-interacting Hamiltonian. It would be inter- 
esting to understand whether these deviations result from our 
limited numerical accuracy in evaluating so small free-energy 
differences, or from additional considerations that affect the 
physics of this system at such large deformations. 

In our elastic Ising model, the low-temperature free en- 
ergy of in-plane positional fluctuations is determined by the 



dimensionless ratio b 



between the magnetoelastic in- 



teraction strength to the lattice rigidity. This ratio sets the 
deformation angle /3 of each isosceles triangle in the ground- 
state (see Eq. ([2]) and Fig. QJ3). In the corresponding system 
of buckled colloids, the deformation angle /3 of each isosceles 
plaquette in the close-packed state (which is equivalent to the 
Ising model's ground state), is similarly dictated by the ra- 
tio of each sphere's diameter to the separation between the 
confining walls [21] ■ We expect the colloidal system to ex- 
hibit a similar order-by-disorder effect that will be governed 
by the deformation angle fJ, and therefore we plot in Fig. [Bp 
the free-energy difference Ao — Ai between straight and bent 
segments vs the geometrical parameter fi rather than vs the 
physical parameter b. 

For extremely rigid lattices that hardly deform (fJ ~ 60°), 
straight and zigzagging stripes are almost equivalent in terms 
of their particle displacements and therefore _4o w A\. As 
f3 increases, straight stripes develop an entropic advantage 
which comes from the fact that the straight-stripe ground 
state is the most anisotropic and hence has the most nonuni- 
form distribution of eigen-frequencies, and thus the maximal 
entropy [13] . For extremely large deformations (/? > 100°) we 
observe a decrease in Ao—Ai which reaches a minimal value at 
fJ ~ 155° and then increases again as /3 — > 180° . In the SI Ap- 
pendix we show how this non-monotonic behavior results from 
the numerical structure of the dispersion relations. In partic- 
ular, we identify in Fig. S10 the region in reciprocal space 
which dominates the free-energy difference between straight 
and bent stripes. The magnitude of the free-energy difference 
in this region grows monotonically with fi, however the size of 
this region decreases, which gives rise to the non-monotonic 
behavior seen in Fig. [Bp for 60° < /3 < 155°. For fJ > 155° 
the free-energy difference has a different wavevector depen- 
dence (see Fig. S10) which gives rise to the second increase 
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in the total free-energy difference at such large deformations. 
It would be interesting to theoretically understand the deeper 
origins and possible implications of this non-monotonic behav- 
ior. Note that for large values of /?, the distance to some of 
the next-nearest neighbors becomes smaller than the distances 
between nearest neighbors. This introduces further complica- 
tions beyond the analysis presented in this paper which as- 
sumes only nearest-neighbor interactions. Overall, our nu- 
merical results show that for all ground-state deformations, 
straight stripes are preferred entropically, however their en- 
tropic advantage is extremely small. We observed the same 
qualitative behavior when allowing the system's total volume 
to vary. 

Before concluding we note that the free energy coefficient 
A given in Eq. Q and plotted in Fig. [6] is defined for each 
specific realization of the stripes. Randomly zigzagging stripes 
which mix straight and bent segments are highly degenerate, 
and their free energy possesses also a configurational entropy 
(or entropy of mixing) which competes with the vibrational 
entropic advantage of straight stripes that we found. However, 
due to the one-dimensional character of the ground state, this 
configurational entropy is sub-extensive in system size and 
scales as y/~N. The vibrational free energy of Eq. enters 
as an extensive quantity which scales linearly with N in T ', 
and is thus dominant in the thermodynamic limit. As a re- 
sult, straight stripes are favored. This may be the reason for 
the limited ability of our simulations to reach perfect straight 
stripes, and for the smaller values of P(iii) / 'P(iv) for smaller 
systems in Fig. [3p. To test this, we repeated the simulations 
at parameter values where the difference in A between straight 
stripes and zigzags is smaller and indeed found a weaker pref- 
erence for straight stripes. 



In summary, our exact microscopic description of this 
previously-studied [5] [B] elastic Ising model reveals a highly- 
degenerate, partially-disordered ground state. We find a much 
richer behavior since entropy lifts the degeneracy at any pos- 
itive temperature. However, large free-energy barriers be- 
tween the ground-state configurations induce a glassy phase 
of zigzagging stripes. Although equilibration is strongly hin- 
dered in this glassy phase, the straight-stripe structures ob- 
tained in it result from the minute entropic differences in the 
equilibrium free energies of the various stripe realizations. The 
fact that this model is amenable to analytic treatment makes 
it appealing as a prototypical model for studying such order- 
by-disorder phenomena in a broader context. Moreover, the 
current approach to experimentally studying frustration re- 
lief by lattice deformations in antiferromagnets is based on 
quite indirect measurements of lattice deformations [451 146] . 
Our work makes direct contact to a colloidal system [22] in 
which frustration and its relief are governed by similar physi- 
cal mechanisms, yet it has the advantage that local deforma- 
tions can be directly measured in it. Interestingly, several re- 
cent experimental works are focused on measuring the normal 
modes of vibration in colloidal systems '4T 48, 49, 50 . On top 
of understanding deformable antiferromagnets and mesoscopic 
model systems for them, we expect that results obtained for 
our system will shed light on questions such as the statisti- 
cal mechanics of sphere packings and the physical origins of 
glassy dynamics. 
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Fig. 1. Model and its ground state: A) Two-dimensional triangular network of antiferromag- 
netic Ising spins connected with harmonic springs. B) Deformation angle (3 vs the dimensionless 
ratio b = between the magnetoelastic interaction strength and the lattice rigidity, from 
Eq. p). Inset: Isosceles triangle with head angle j3, satisfied bonds (thin lines) compressed by 
s and frustrated bond (thick line) stretched by /. C) The five possible configurations (up to 
rotations and spin inversions) of a particle and its neighbors, with a single frustrated bond in 
each triangle. 




Fig. 2. Zigzagging stripes: Ground-state configurations are generated by stacking layers of 
particles with alternating spin and with arbitrary relative polarities between successive layers. 
The figure shows only simple configurations with periodic sequences of the straight (Hi) and 
bent (iv) segments, as defined in Fig. However ground-state configurations do not neces- 
sarily have a finite unit cell. A) Straight stripes, for which all particles are in the state (Hi). B) 
Bent stripes, for which all particles are in state (iv). C) Zigzagging stripes with a more com- 
plicated unit cell comprised of both (Hi) and (iv). The shaded regions represent the unit cells 
used for the low-temperature expansion explained in the text. (A), (B), and (C) have one, two 
and five particles per unit cell, respectively. Thick blue lines represents the stretched frustrated 
bonds and thin blue lines the compressed bonds for which the antiferromagnetic interaction is 
satisfied. 
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Fig. 3. Monte-Carlo simulations: A) Energy E per particle above the ground-state energy 
Eqs vs temperature T for different cooling rates r as indicated in the legend. System size 
is TV = 6400. Similar results were obtained for TV = 1600 and TV = 25600. B,C) 
Probabilities of finding the local configurations (Hi) and (iv) defined in Fig. at T = 
following cooling at different rates V for various system sizes TV, as indicated in the legend. 
Note the logarithmic scale for the ratio P(iii) / P(iv) in (C). Error bars are smaller than the 
symbols. Model parameters in all simulations are J = 1, e = 2, K = 8, (X = 1. Thus 
b = -j^- = 0.25. This yields a deformation angle of /3 = 86° in the ground state and a 
difference in the free-energy coefficient, Eq. J4), of f/^40.05 between straight and bent stripes. 
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Fig. 4. Ordering with decreasing cooling rate: Portion of the system in its final configuration 
after cooling at rates r = 10~ 4 (A), 10~ 6 (B), 10~ 8 (C). Model parameters are J = 1, 
e = 2, K = 8, O, = 1. Thus, b = = 0.25, and in the ground state each triangular 
plaquette is deformed to an isosceles with head angle j3 = 86°. System size is TV = 25600. 
For the fastest cooling rate, the system falls into a disordered state (A), for which the fraction of 
particles in the local configuration of a straight stripe (see Fig. UXD is P(iii) = 0.31 and the 
fraction in that of bends is P(iv) = 0.52. For the slower cooling rates, zigzagging stripes are 
formed (B,C). The blue lines are guides to the eye, which highlight a line of frustrated bonds, 
and emphasize that as the cooling rate is decreased from (B) to (C), the stripes become more 
straight. In (B), P(iii) = 0.54 and P(iv) = 0.44, whereas in (C), P(iii) = 0.9 and 

P(iv) = 0.1. 
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Fig. 5. Examples for some of the unit cells used in evaluating the free energy in Fig. \§\ All 
unit cells with p s = 1/3 containing 3 (A), 6 (B), and 9 (C) particles, and with p s = 1/2 
and 4 (D) and 8 (E) particles. 
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Fig. 6. Normalized entropic contribution, Eq. @ vs fraction p s = P (Hi) / [P (Hi) + 
P(iv)] of straight-stripe segments for deformation angles (3 = 70° (A), 100° (B), and 
155° (C). A is normalized by its extreme values for straight stripes A,\ = A.(ps — 1) 
and for bent stripes A.Q = >A\p& = 0). The gray lines indicate the prediction of a one- 
dimensional effective Hamiltonian of non-interacting straight and bent segments. The size of 
the unit cell corresponding to each ground-state configuration is indicated in the legend. We 
calculated the free energy for all the 95 distinct unit cells consisting of up to 10 particles. Except 
for straight (p a = 1) and bent (p s = 0) stripes, these have ^ < p s < ^. Additionally, 
of the 363 unit cells with 11-13 particles, we calculated the free energy for 19 unit cells with 
extreme values of ps- In numerically evaluating the integration over the Brillouin zone, we 
require a relative numerical accuracy smaller than 10 — ^ , and therefore for (3 = 155° omitted 
from the plot 19 unit cells (with 11-13 particles) and for /3 = 100° omitted one unit cell (with 
13 particles). See SI Appendix for more details on the calculations and the results. D) Entropic 
advantage of straight stripes (p s = 1) over zigzags (p s = 0) vs deformation angle. The blue 
circles indicate the angles for which results for larger unit cells are given in (A-C). 
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